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BaBi03 is a prototype "charge ordering system" forming 
interpenetrating sublattices with nominal valence Bi 3+ and 
Bi 5+ . It can also be regarded as a three-dimensional version 
of a Peierls insulator, the insulating gap being a consequence 
of an ordered distortion of oxygen atoms. When holes are 
added to BaBiC>3 by doping, it remains insulating until a 
very large hole concentration is reached, at which point it 
becomes superconducting. The mechanism for insulating be- 
havior of more lightly-doped samples is formation of small 
polarons or bipolarons. These are self-organized point defects 
in the Peierls order parameter, which trap carriers in bound 
states inside the Peierls gap. We calculate properties of the 
polarons and bipolarons using the Rice-Sneddon model. Bipo- 
larons are the stable defect; the missing pair of electrons come 
from an empty midgap state built from the lower Peierls band. 
Each bipolaron distortion also pulls down six localized states 
below the bottom of the unoccupied upper Peierls band. The 
activation energy for bipolaron hopping is estimated. 
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I. THE MODEL 

Pure BaBiC>3 has a distorted perovskite structure jt). 
The nominal valence of Bi is 4+, leaving a single s- 
electron per Bi atom. If BaBiC>3 had perfect cubic per- 
ovskite symmetry, then (in independent electron picture, 
or band theory) it would be metallic, with the Fermi level 
half into a broad band made of Bi 6s electrons (partially 
antibonded with O 2p a states.) The large size of the Bi 
6s orbital indicates that single-electron (band) approxi- 
mation should be good 0. Therefore the origin of in- 
sulating behavior should lie in the structural distortions 
which double the unit cell. The simplest interpretation 
is that BaBiC>3 is a prototype d=3 Peierls insulator with 
a simple lattice dimerization. 

BaBiC>3 can be hole-doped by Pb substitution for Bi 
or K substitution for Ba. At a quite high critical concen- 
tration x c of holes per cell (x c ss 0.65 for Pb or w 0.35 
for K), an insulator to metal transition occurs, with the 
metal a superconductor. The high superconducting tran- 
sition temperature (T c < 30K) is compatible with a con- 
ventional electron-phonon mechanism, but neither theory 
H nor spectroscopy [Q is able to make a convincing con- 
firmation. Electron-phonon effects certainly play an im- 
portant role in the insulating part of the phase diagram. 



This paper explores the self-trapped polarons which oc- 
cur at low doping and which provide the most plausible 
mechanism explaining how these materials remain insu- 
lating to such high hole concentrations x c . 

We use a simple and yet fairly realistic model which 
gives a microscopic description of the Peierls distortion. 
The model is illustrated in Fig. |l]. There is one elec- 
tronic degree of freedom per Bi atom, the amplitude of 
the 6s electron orbital at that site, with corresponding 
creation operator c e . The index I is a composite, standing 
for {I, a) where the integer vector I = (l Xl l y ,i^) locates 
Bi atoms on a cubic lattice of lattice constant a=4.28A, 
and cr = ±1 is the spin index. The Bragg vector of 
the dimerized lattice is Q = (tt,tt,tt). The phase factor 
exp(iQ • t) = (— 1) L , where L = £ x + £ y + £ z , separates 
the simple cubic lattice into two sublattices, A and B, 
with (— 1) L equal to 1 on A and -1 on B. The 6 nearest 
neighbors of each atom lie on the opposite sublattice. 
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FIG. 1. Two-dimensional section of the Rice-Sneddon 
model. The small filled circles are oxygen atoms which are 
displaced by amounts u from the midpoints of the bonds. 
The central Bi atom (labeled /) is on the A sublattice (nomi- 
nally Bi J+ ) where electrons accumulate because oxygens have 
moved away, lowering the potential energy of electrons. 

In each cell there are three vibrational degrees of free- 
dom, the oxygen displacements along the Bi-0 bond di- 
rections. Vibrations perpendicular to bond directions are 
assumed not to influence electrons. The oxygen displace- 



1 



ments u{£, a) are labelled by indices a — x,y,z which 
refer to the oxygen located at position a£ + (a/2)a. It is 
convenient to define a local dilation or breathing ampli- 
tude e(£) on the £-th Bi atom, 

e{?)= H£,+,a)-u(£,-, a)] (I) 

a— x,y,z 

where u{£ , ±, a) is the displacement of the oxygen atom 
at a£ ± (a/2)d. Thus u(£, +,a) is the same as u(£, a), 
while u{£, — , a) is the same as u(£ — d, a), that is, it is 
the displacement of the oxygen atom located at — (a/2)d 
relative to the Bi atom I. For a positive distortion e(£), 
the Bi 6s electron has its energy lowered by —ge{£) where 
g is the deformation potential. Each oxygen displace- 
ment costs elastic energy (K/2)u(£,a) 2 . Using the Ra- 
man measurement of fko=70 meV for the oxygen "breath- 
ing mode" 0, the value of K is 19 eV/A 2 . The resulting 
Hamiltonian is 

H = - t E c k'-5E e W c ^ + ^E u ^ a ) 2 

<l,Z>> i i Q 

(2) 

where the hopping summation index <£,£'> goes over 
each nearest neighbor pair both forward and backward, 
and conserves spin. The variables e(£) for different £ are 
not independent, only the original variables u{£, ±, a) of 
Eq. (|l|) are independent. Piekarz and Konior |?]] have 
introduced new decoupled breathing mode variables, at 
the expense of more distant terms in the coupling (second 
term) of Eq. (§). 

For the time being we make the adiabatic approxi- 
mation, treating oxygen mass M as infinite (neglecting 
oxygen kinetic energy.) The electronic band structure is 
modelled with a simple (ssa) hopping integral t between 
nearest neighbor Bi atoms. When there is no oxygen 
displacement, this yields the dispersion 

e(fc) = —2t(cos(k x ) + cos(ky) + cos(fc 2 )) (3) 

Comparison with band theory |3|J^-[ll|] indicates that the 
value t = 0.35 ± 0.5 eV is appropiate for the band-width 
12t w 4 eV of the conduction band in cubic BaBiC^. Out- 
values for the parameters are summarized in table S. 



TABLE I. Estimated values of parameters for BaBi03 . All 
parameters have «10% uncertainty. 





Fitted Parameters 




t 




0.35 eV 


K 




19 eV/A 2 


9 




1.39 eV/A 




Derived Parameters 




r = g 2 /Kt 




0.30 


S = Tiuj/t 




0.017 



If energies are measu red i n units of t and lattice dis- 
placements in units y^t/K, then the model contains 
only a single dimensionless coupling constant yT 7 where 
r = g 2 /Kt. The dimensionless Hamiltonian H/t is 
just Eq.(|^) with the substitutions t — > 1, K — > 1, and 
g — * VT. Later when non-adiabatic effects are treated, a 
new dimensionless energy ratio 6 = hto/t enters. 

Because hopping goes only between inequivalent sub- 
lattices, the band energy has the symmetry property 
e{k + Q) = — e{k). The bands are symmetric around 
energy 0, which is the Fermi energy at half filling. There 
is perfect "nesting" at this value of the Fermi energy, 
because for every state with ek — 0, there is another 
state at wavevector k + Q which is also at the Fermi 
energy. This perfect nesting property is an artifact of 
the nearest neighbor hopping, and is destroyed by next- 
neighbor hopping terms. Because of perfect nesting, it 
is particularly simple to find the ground state of the 
undoped (half-filled) case, obtaining a BCS-like integral 
equation for the Peierls gap 2 A, which has a solution 
A « 6.5texp(-0.291/r) for weak coupling (r < 1). 
This was discussed in a previous paper |l^|. The next 
section will treat this problem in the extreme atomic 
limit (r> 1) where a particularly simple approximation 
A = 12Ft applies. 

An important consistency check on the Hamiltonian is 
the fact that three parameters (t, K, and g) fit four inde- 
pendently determined properties of BaBiC>3. The value 
5=1.39 eV/A is chosen to yield a Peierls band gap 2 A = 
2 eV, agreeing with optical measurements |l3| . Then the 
oxygen sublattice displacement uo — A/65 is predicted 
by the model to be 0.12 A, while experiment sees 0.09 A 
§0- 

Our model Hamiltonian contains a good explanation 
of the difficulty of making free carriers: doped-in holes 
self-localize, making small polarons and bipolarons. In a 
future paper we will re-examine the case of heavy dop- 
ing, previously treated by Yu, Chen, and Su pa]. Here 
we concentrate on dilute doping. We have previously [T^ ] 
looked at the criteria for polaron formation in the case of 
dilute doping into an empty band (no Peierls distortion) , 
as occurs in materials like BaSnOs. The critical coupling 
strength T C (P) for polaron formation was found to be 
1.96, and F C (B) for bipolaron formation was found to be 
0.99. We have also studied the case of dilute doping into 
the half-filled band by an analytic variational approxima- 
tion [ fn| . Critical coupling strengths were dramatically 
reduced, to T C (P) = 0.18 and T C (B) = 0.15. Here we 
provide detailed numerical pictures of the behavior of 
lightly doped BaBiOa, including predictions about the 
optical spectrum. In a companion paper \ \n\ we examine 
the self-trapped exciton state which forms when undoped 
BaBiOs is excited optically. 
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II. EXTREME ATOMIC LIMIT T > 1 

The value T ~ 0.30 is appropriate for BaBi03. This is 
intermediate between weak and strong coupling. Beyond 
its applicability to BaBiC>3, the model is intellectually 
interesting in its own right, and we can be curious how its 
behavior evolves with coupling constant T over the whole 
range to oo. We now solve the problem in the extreme 
atomic limit T — > oo which just means that we ignore 
the hopping term (f — » 0), simplifying the mathematics 
greatly. The results yield insight into the intermediate 
region T ~ 0.30. 



A. Bipolaron crystal 

When the Bi s level is half-filled, the strong-coupling 
solution consists of putting two electrons on each Bi atom 
on the A sublattice and no electrons on each Bi atom on 
the B sublattice. This is the extreme atomic limit of a 
Peierls charge density wave. The paired electrons can be 
given various names: "bipolarons" or "negative U cen- 
ters." The mechanism stabilizing these objects is oxygen 
motion. Each oxygen has as nearest neighbors one A 
and one B sublattice Bi atom. Therefore it experiences 
a force away from the more negative A site and toward 
the more positive B site. The piece of the Hamiltonian 
(0) containing this force is 



gu{i, a 



(4) 



where £' — £ + a is the other member of the pair of near- 
est Bi atoms. Let us define the Peierls order parameter 
(amplitude of charge density wave) po by 



£<c|c £ >=l + (-l) J 



Po- 



(5) 



In the undistorted cubic structure, po = 0. In strong 
coupling, po — > 1, that is, the Bi atoms have charge 
2 on sublattice A and on sublattice B. The force is 
2gp on each oxygen, directed toward its nearest B-type 
Bi neighbor, counteracted by a harmonic restoring force 
Ku, yielding an optimum displacement uq — 2gpo/K 
or a "breathing amplitude" e(£) = (— l)Xo, where, the 
oxygen breathing order parameter eo is V2gpo/K. These 
relations remain true for arbitrary coupling (0 < pa < 1). 
In strong coupling, site energies of Bi s orbitals are thus 
— 12g 2 /K on A sites and +12g 2 /K on B sites. For hop- 
ping t = 0, the lowest energy excitation is to move one 
electron from an A site to a B site. This costs energy 
2A = 2Ag 2 /K = 2ATt which is the value of the Peierls 
gap. We shall use the notation A to denote the value 
of the gap in the atomic limit, and A for the value for 
the actual coupling strength under consideration. The 
limiting answer 2Ao = 24g 2 /K can also be found from 
the integral equation (23) of ref. [O, by taking A » t. 



We get a measure of the value T c where the crossover 
occurs between strong and weak coupling by comparing 
total energies. In the atomic limit, this is —6g 2 /K per 
atom (— 12g 2 /K per electron of on-site orbital energy, 
but at a cost of +6g 2 /K per atom of elastic energy.) In 
the opposite (<?=0) limit, the electrons gain on average 
— 2t of derealization energy by band formation in the 
undistorted crystal. These energies are equal at T c = 
1/3. 

So far we are safe in the adiabatic approximation, since 
the vibrational energy huj is small compared with the 
Peierls gap in strong coupling. However, the true lowest 
energy excitation will cost less than 2A once we allow 
relaxation of the atoms. Thus the optical spectrum will 
have non-adiabatic Franck-Condon character, with min- 
imum excitation energy corresponding to a "self-trapped 
exciton." The companion paper |l7f pursues this topic. 
The rest of this paper is devoted to the case where one 
or two holes are doped into the half-filled band. We will 
consider both the ground state and the new excitations 
induced by the doping. 
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FIG. 2. Schematic structures: (a) is the undistorted high 
temperature structure; (b) is the low temperature Peierls 
phase of pure BaBiOa; (c) has a hole polaron in the central 
cell; and (d) has a hole bipolaron. Circles denote Bi atoms, 
with filled circles being Bi 3 + ions with two s electrons, shaded 
circles being Bi 4+ ions with one s electron, and open circles 
being Bi 5+ ions with no s electrons. X's denote oxygen ions. 



B. One hole in the atomic limit 

To remove an electron, keeping the Peierls distortion 
fixed, costs energy 12g 2 /K (the amount of Peierls at- 
tractive energy lost.) When t — 0, the hole sits on 
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some particular A sublattice site as a singly charged Bi 4+ 
ion. The 6 surrounding oxygens now experience a dimin- 
ished repulsive force, g instead of 2g, toward their Bi 5+ 
neighbors. The local Peierls displacement u\ is there- 
fore g/K, or half of uq. This reduces the elastic energy 
by 6 times K(uq — u\)/2 or 9g 2 /K but does this at the 
expense of raising the energy of the 6s orbital at this 
site from —12g 2 /K to —6g 2 /K. Since this orbital is oc- 
cupied once, the net energy saved by this local distor- 
tion is 3g 2 /K. The resulting object, a point defect in 
the Peierls order parameter, is a small hole-polaron. In 
a system with N Bi atoms, the isolated polaron corre- 
sponds to the system having N — 1 electrons. Its energy 
E[P] = E gB (N - 1) - E gs (N) is 9g 2 /K, smaller than 
the removal energy 12g 2 /K when no relaxation is per- 
mitted. The difference is the polaron trapping energy, 

/••V/' ■•'-'/•; K. 

Photoemission would measure the polaron formation 
energy 9g 2 /K rather than \2g 2 /K provided the removal 
were done slowly enough that oxygen atoms could move 
to their new optimum positions. In fact, this is surely not 
the case, and one expects to see instead a peak photoe- 
mission intensity at the adiabatic energy 12g 2 /K , with 
a series of Franck-Condon vibrational sidebands with 
a Gaussian envelope extending down to the threshold 
9g 2 /K. A detailed discussion of the corresponding ef- 
fect in LaMnC>3 was given earlier |L8| . 

The hole polaron is shown schematically in part (c) of 
Fig. H It is clear that the orbitals on the surrounding 6 
Bi 5+ ions must also be altered by the oxygen relaxation. 
Rather than having their energy pushed up by 6guo, their 
energy is pushed up by 5guo + gu\, or \\g 2 / K rather 
than 12g 2 /K. Thus the polaron defect has an interesting 
spectrum of gap states. There is the singly occupied A- 
sublattice state located at A/2 above the lower Peierls 
band at — A, and there are six empty -B-sublattice states 
located at energy A/ 12 below the upper Peierls band at 
A. This spectrum is shown on the left side of Fig. ||. 
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bipolaron 



FIG. 3. Energy spectrum of the Rice-Sneddon model in the 
atomic limit (t — 0) with one hole present (left side labeled 
"polaron" ) or with two holes present (right side labeled "bipo- 
laron.") In both cases, the valence band contains N/2 — 1 
doubly-filled states at energy —A, and the conduction band 
has N/2 — 6 degenerate states at energy A, with 6 additional 
empty states bound below them. 



C. Two holes in the atomic limit 

Two polarons spatially separated cost two times 
9g 2 /K, but if they occupy the same local site (with op- 
posite spin, 5 = 0) then it is easy to see that the cost 
is only 12g 2 /K. This is because removal of two elec- 
trons costs 2Ag 2 /K in lost Peierls energy, but when sit- 
ting on a single A site (converting the Bi 3+ to a Bi 5+ ), 
the local strains on the six surrounding oxygens all dis- 
appear (u 2 — 0.) This gains back \2g 2 /K that was 
previously paid in strain energy. The local s orbital 
at the site of the bipolaron now has energy because 
of no local strain, but unlike the polaron case, the loss 
of Peierls energy of this state does not alter anything 
since the level is now empty. Thus the total energy 
of the bipolaron is E[B] = \2g 2 /K. Relative to free 
holes which cost 24g 2 /K, the bipolaron trapping energy 
is E t [B] = \2g 2 /K. The energy released when two po- 
larons bind into a singlet bipolaron is E\, = 6g 2 /K. How- 
ever, this neglects on-site Coulomb repulsion U. We 
expect a reduction of the bipolaron binding energy to 
Eb = Qg 2 /K — U. The Hubbard U is defined as the en- 
ergy difference between charged ions E(Bi 3+ ) + _E(Bi 5+ ) 
and "neutrals" 2E{B\ i+ ) for BaBiC>3 with oxygen atoms 
frozen in cubic perovskite positions. Experimentally, the 
absence of free spins in doped BaBiOs indicates that Bi 4+ 
(with S — 1/2) is disfavored. Therefore, U cannot be 
large enough to destabilize bipolarons. Vielsack and We- 
ber Q have done constrained density functional calcula- 
tions to obtain U. They fit their results to an effective 
U of 0.6 ± 0.4 eV of a single band model like ours. This 
is similar in size to our attractive energy of 6g 2 / K i=s 0.6 
eV. When hopping is turned back on, the attractive en- 
ergy goes down, but the bipolaron also spreads out, so 
that the repulsive U is less effective. Thus U is small 
enough to permit bipolarons to exist. 

The spectrum of BaBi03 with a bipolaron consists of a 
filled lower Peierls band containing N/2 — 1 states at en- 
ergy —A = —12g 2 /K, one empty midband state E = 0, 
and N/2 — 6 upper Peierls band states at energy +A. Six 
new empty states appear below the empty upper Peierls 
band, at energy 5A/6. These are the s orbitals on the six 
B sites surrounding the bipolaron, which have only 5 out 
of 6 first neighbor oxygens with displacement uq, and one 
undisplaced. The spectrum is shown on the right side of 
Fig. | 
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D. corrections for small hopping t 

For the perfectly ordered half-filled band, the Peierls 
gap diminishes quadratically as t — Ao/12r is turned on, 

A = Ao(l-^ + .-). (6) 

This follows by expanding the integral equation (23) of 
ref. for small values of t/ A. The amplitude p of the 
CDW diminishes the same way, from 1 to 1 — l/48r 2 + - ■ ■. 

The eigenstates are now ex tended Bloch states of en- 
ergy X(k) = ±-\/A 2 + e(k) 2 , where e(k) is the eigen- 
value of the undistorted lattice, — 2t ^ Q cos(k a ). Thus 
the least expensive delocalized hole states are at the 
top of the lower Peierls band where e(k) — and cost 
(12g 2 /K)(l - l/48r 2 ). The competing bipolaron state 
can no longer be completely localized on a site because 
the charge on the nominal Bi 3+ ions is 1 + p < 2. It is not 
easy to develop a systematic 1 /Y expansion for the local- 
ized solutions. The problem is that as oxygens move and 
the polaron localizes, wavefunctions and charges of all oc- 
cupied delocalized states must be altered self-consistently 
near the polaron. 

III. NUMERICAL STUDY OF POLARON 

The algorithms used are exp lained in detail in the MA 
thesis of the first author |l9[ . A large supercell (typi- 
cally between 100 and 1000 Bi atoms) was chosen, with 
periodic boundary conditions applied. It was useful to 
have the translation vectors A, B, C of the supercell not 
lie along symmetry axes, but to use small asymmetrical 
deviations in order to lift degeneracies. The translation 
vectors all translate A sublattice points to other A sub- 
lattice points, so that the Peierls state fits the cell. Trial 
oxygen positions are chosen, the Hamiltonian matrix is 
diagonalized, and forces on oxygen atoms are calculated 
by summing over occupied wavefunctions to get the elec- 
tronic charges on the two adjacent Bi atoms. A vari- 
able metric method was used to iterate towards op- 
timal oxygen displacements. We confirmed that for the 
half-filled case, the density of states, oxygen coordinates, 
and Peierls gap all agreed well with calculations by other 
methods. The half-filled ground state serves as the "vac- 
uum" for the rest of our work. To study doping by one 
or two holes per supercell, a cell with «200 atoms is suf- 
ficient if the coupling constant Y > 0.22. For smaller Y, 
larger cells are needed. 




FIG. 4. Trapping energy of a polaron (upper panel) and 
a bipolaron (lower panel) computed numerically. Bold solid 
lines are strong-coupling predictions. Thin lines are approxi- 
mate 1/F 2 corrections discussed in Sec. III. A. 

A. total energy 

By simply restricting the occupancy to be N — 1 or 
N — 2, we obtain solutions with polaron or bipolaron lo- 
cal defects. The numerical results for the "trapping ener- 
gies" are shown in Fig. [l| Stable polaron solutions occur 
whenT > Y C (P) as 0.18. Stable bipolaron solutions occur 
when r > Y C (B) ss 0.15. At lower values of Y an inserted 
hole just stays at the top of the valence band, causing an 
order 1/N renormalization of the Peierls gap. For T > r c , 
such states were always found to be numerically unstable 
toward formation of local defects. The onset of localized 
solutions seems to occur continuously, unlike the case 
where an empty s-band is doped and localized solutions 
appear discontinuously (in adiabatic approximation) at 
much larger values Y c — 1.96. There are several reasons 
for the ease of polaron formation in the Peierls state rela- 
tive to the empty band. Perhaps foremost is the fact that 
there is less to be gained by derealization, because low- 
energy delocalized states are occupied. Also, the Peierls 
state has lattice strain already built in, so polaron forma- 
tion can occur by reduction of the Peierls lattice strain 
rather than by initiation of new lattice strain. 

Although we do not have a systematic 1/r perturba- 
tion theory, nevertheless, a very simple correction to the 
r — > oo expression can be made. The trapping energy 
of polarons E t [P] is the difference between the energy A 
to create a hole without relaxation, and the energy E[P] 
of the relaxed polaron. In the r — > oo limit, this gives 
E t [P]/t = 12Y - 9Y. The first term, A /t, the energy to 
put a hole into the "vacuum" without lattice relaxation, 
can be corrected by Eq. (|) to A(Y)/t w 12r~ 1/4I\ We 
do not know the form of the correction to ], but find 
that the correction is weaker than the vacuum correction. 
Therefore, correcting only the vacuum term gives the for- 
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mula E t [P]/t « 3r~ 1/4T. The r -> oo formula (labeled 
t — > 0) and the vacuum-corrected formula are shown in 
Fig. H as thick and thin solid lines. The fit to the po- 
laron is very good down to T = 0.5. For the bipolaron, 
the corresponding formulas are E t [B] = 2 A — E[B], and 
E t [P]/t w 12r - 1/2Y. The vacuum-corrected fit is not 
as good as for the polaron. 



B. size and shape of polarons 

In the t — > limit, an inserted hole sits on a single site. 
Hopping causes the hole to spread, gaining delocalization 
energy. The amount of positive hole charge p(f) at a site 
r is defined as 



OCC,P 



p(f)= J2 \^^c(r)\ 2 -J2 \*iA?)f 



(7) 



where "vac" and P indicate the single particle states ^> 
for vacuum and polaron cases. In Fig. [5], the hole charge 
at the central site r — is plotted. For T > 0.5, the 
polaron is almost completely on the central site, whereas 
for the bipolaron, the available charge on the central A- 
type atom is essentially all depleted for T > 0.4, but 
until a somewhat larger T, some charge must be removed 
from surrounding sites. At the realistic value T = 0.3, 
the polaron depletes 0.84 electrons from the central A 
site, and the bipolaron depletes 1.74 electrons, which is 
almost all of the 1.80 electrons which occupy A sites in 
the vacuum state. 
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FIG. 5. The total charge removed from the central site of 
the hole polaron and bipolaron. The solid line is the total 
charge at that site in the Peierls state, the maximum available 
for removal. 



gain electrons when a polaron is formed nearby; neigh- 
boring A sites are more depleted of electrons than if B 
sites had been unaffected. The distance-dependence of 
the hole charge is shown in Fig. (|^) to be exponential, 
with decay length 0.39a. The corresponding decay length 
for the bipolaron disturbance is 0.32a. 
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FIG. 6. Charge versus distance for the hole polaron when 
T = 0.30. The calculation used a cell with AT=1000. Sites 
with depletion of electron density (charge Sp > 0) are denoted 
by circles; sites with excess electron density (Sp < 0) are 
denoted by squares. Beyond r = 5a, overlapping charge from 
neighboring cells affects the answer. 

If small polaron formation occurs in continuous fashion 
for r > r c , as suggested by Fig. |5| shows the continu- 
ous evolution with r, for T > T c , from large CDW-like 
(bi) polarons to small (bi) polarons. The decay length A 
(defined as the negative reciprocal of d\n(8p)/dr) should 
diverge as T is reduced to T c . This cannot be verified 
numerically because of the small linear size of available 
computational cells in 3 dimensions. We do see a trend 
in the right direction, shown in Table O. 

TABLE II. Decay constants for the hole charges for a 
(bi)polaron 
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1.97 


0.6 


0.24 


0.98 


0.23 
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0.84 
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The polaron defect diminishes the amplitude of the 
Peierls charge-density wave in its local vicinity. This is 
seen by examining the charge on atoms near the polaron. 
At r = 0.3, to make a polaron, 0.16 electrons are taken 
from sites other than the central site. We find that neigh- 
boring B sites do not lose electron density, but rather 
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C. Conduction band edge states 

The polaron defect perturbs the band Pcicrls states in 
a local region. This perturbation is sufficiently strong 
to create new localized states at the bottom of the con- 
duction band. This effect is particularly simple in the 
strong-coupling limit, where we already discussed (see 
Fig. and nearby text) six states bound by A/f 2 (near 
a polaron) or A/6 (near a bipolaron). These states can 
be classified as one s-like (Ai g ), a c?-like doublet (E g ) and 
a p-like triplet (T± u ). Therefore, as hopping t is turned 
on, we can expect the six bound states to split into sin- 
glet, doublet, and triplet, and for weak coupling, perhaps 
to become unbound. The numerical results are shown in 
Fig. 0. The s-like state is bound only for T > T c (s) w 
0.30 (for a polaron) or 0.28 (for a bipolaron), while the 
other states are nearly degenerate and look to persist all 
the way to the critical T for polaron or bipolaron forma- 
tion. Numerical finite-size effects prevent high accuracy 
in these estimates. The reason for the weaker binding 
of the s-like state is that it couples to the polaron state 
which splits off from the valence band, forming a spread 
out polaron with no radial node, and an empty state 
bound below the conduction band with a radial node. 



nately, for exciting the mid-gap bound hole to the valence 
band) is zero in first order at the k = (tt/2,tt/2,tt/2) 
band edge. There is a non-zero matrix element to the p- 
like states split off from the bottom of the upper Peierls 
band, with energy w 1.88 A or 1.88 eV. These numbers 
should not be taken as reliable predictions of actual spec- 
tral lines for several reasons. First, the predictions are 
sure to be somewhat sensitive to the model, especially 
to the fact that our model omits the disorder potential 
caused by dopant atoms. Second, there will be significant 
Franck-Condon effects coming from the fact that if the 
empty midgap states are suddenly occupied by an optical 
transition, local oxygen relaxations can lower the energy 
of these states. Finally, as explained in the companion 
paper, the Peierls gap of pure BaBiC>3 in this model is 
not "clean." When an electron-hole pair is created in 
the pure material, self-trapping of the electron and hole 
should occur, causing absorption in Franck-Condon side- 
bands below the energy of the nominal Peierls gap 2 A. 
This may obscure additional spectral features introduced 
by bipolarons formed by doping. 

IV. MOBILITY OF BIPOLARONS 



<l 

Qj 




FIG. 7. Energies of localized single particle states lying in 
the gap. The lower panel shows the states which derive from 
the valence band; the polaron state is occupied once (spin 
degeneracy = 2.) All other states are empty in the ground 
state. The upper panel, showing states localized below the 
conduction band edge, has a different scale, both horizontally 
and vertically. 

In principle, these bound states can be seen spectro- 
scopically, as weak mid-gap absorption proportional to 
doping concentration. We believe that bipolarons are the 
stable point defects. At r=0.3, the doped-in holes are 
bound at s» 0.7A or 0.7 eV above the top of the lower 
Peierls band. The optical matrix element for exciting 
an electron from the valence band into this state (alter- 



At T=0 K, band formation occurs, in principle, with 
band width exponentially reduced by Huang-Rhys fac- 
tors 21.23] exp(— S) coming from vibrational overlap in- 
tegrals. Let us make a crude estimate for the polaron. 
To simplify, we work in the strong coupling limit where 
the polaron charge is all at one site, £. The state \£P > 
representing the polaron at this site, is a product of a 
hole wavef unction eg | vac > (where |vac > is the elec- 
tronic vacuum) times a vibrational wavefunction \e±(£) > 
(the vibrational ground state with six oxygens in altered 
breathing positions around site £.) We want the effec- 
tive hopping matrix element from this state to a degen- 
erate polaron state \£'P > located at a neighbor site £' . 
The Hamiltonian matrix element < £P\TL\£'P > factor- 
izes into 



teS —< vac|c|Hcf|vac >< ei(£)\ei(£ ) 



> 



(8) 



The first factor of Eq. (j|) is the matrix element to 
hop from £ to both being on the A sublattice (since 
that is where the polarons reside.) This is a second- 
neighbor hop on the underlying simple cubic lattice. Our 
Hamiltonian omitted second-neighbor hops, but can now 
be extended to include such a hopping matrix element t', 
smaller than the first-neighbor hopping t. Alternately, 
we could do second-order degenerate perturbation theory, 
which would be rather messy and would probably yield 
something of the same magnitude. 

The second factor of Eq. (j^) is the "Huang-Rhys fac- 
tor," namely the overlap between vibrational ground 
states of two different lattice configurations. Twelve oxy- 
gen atoms change locations between the two states. Six 
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oxygens around the original polaron site change positions 
from u\ to wo, and six oxygens around the new polaron 
site change positions from uo to u\. The corresponding 
overlap parameter is 



exp (-S) = exp[-12(M(j/h)(u - "i)7 4 ] 



(9) 



The exponent, rewritten as S = (A/7iu;)(A/t)(l/48r), is 
approximately 2.8. The polaron band width is thus of 
order 0.06 eV, if we use the band-width 4 eV, and as- 
sume that the second-neighbor hopping is approximately 
t' = t/4. The disorder caused by the dopant atoms is un- 
doubtedly larger than this, so we should expect polarons 
to be Anderson-localized. 

The bipolaron band-width will be much narrower. For 
two electrons to hop to second neighbors is a higher or- 
der process, and the Huang-Rhys exponent S is larger 
by 4 since oxygens move the full amount uq rather than 
by \uq — u\\ = uq/2. Thus the Huang-Rhys factor for 
bipolarons is smaller by 2 x 10~ 4 than for polarons, and 
we expect bipolarons to be completely immobile at T=0 
K. 

At high temperature (fc^T sa TiuS), bipolaron mobility 
should reduce to the classical version, attempt frequency 
times thermal activation factor. Assuming that thermal 
activation is given by exp(— Es/ksT), with E$ the en- 
ergy of the saddle point, we can make a plausible argu- 
ment for how mobility should behave at high T. The sad- 
dle point is illustrated in Fig. |[ In this state, the bipo- 
laron has split into two polarons located at A-sublattice 
first neighbors, with full lattice relaxation to minimize 
the energy at the saddle. This state lies above the bipo- 
laron state in energy by the energy required to unbind 
a bipolaron into a first-neighbor polaron pair, A/2 « 
0.5 eV in strong-coupling approximation. For the ac- 
tual value r=0.3, we estimate the activation energy to 
be Eg ~ A/3 w 0.3 eV. Experimental measurements J2^] 
vary from 0.17 eV to 0.27 eV. 






FIG. 8. Hopping of a small bipolaron through an interme- 
diate state with two polarons. Filled circles are Bi 3+ ions with 
two electrons, and open circles are Bi 5+ ions with no electrons. 
The left panel has a hole bipolaron in the small square, which 
evolves into two hole polarons in the middle panel (shaded 
circles, medium size squares.) One of the polarons then col- 
lapses back to a Bi 3+ ion (right panel, large square) and the 
bipolaron has moved to a new location (small square). 
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